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Abstract 

. An intrinsic metric tensor, a flat connexion and the correspond- 

0^ . ing distance-like function are constructed in the configuration space 

^ . formed by the velocity field and the thermodynamic variables of an 

D . inviscid fluid. The kinetic-energy norm is obtained as a limiting case; 



all physical quantities are Galilean invariant. Explicit expressions are 
given for the case of an ideal gas. The flat connexion is not metric- 
compatible. These results are achieved by applying the formalism of 
• statistical manifolds [Q, ^ to the statistical mechanics of a moving 

^ : fluid, 
o ■ 

CN ' Submitted for publication. 

1 Introduction 

O \ Typically, the configuration of a simple fluid at a certain time t is described by its 

^ . density field p{x,t), its temperature field T{x,t) and its velocity field u{x,t) where x 

^ is the position vector in a (i-dimensional space. This (i- dimensional space is often the 

^, everyday three-dimensional Euclidean space or a two-dimensional surface embedded 

' in three-dimensional space, e.g., a plane or a sphere. In that dimensional space, it 

is possible to calculate the square of the distance between two nearby positions, say 
x and X + dx, in the form of an expression 

Mx^{x)dx^dx'' (1) 

where M\^{x) is the metric tensor, dx^ is the A-component of dx in some coordinate 
system, 1 < A and u < d, and the convention of summation over repeated lower and 
upper indices holds. For example, on a two-dimensional {d = 2) sphere of radius r, 
in the coordinate system defined by the angles and (p (latitude and longitude), the 
non-vanishing components of the metric tensor are 

M^^ = r\ (2) 
M^^ = r^smH, (3) 
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and the corresponding volume element is y det {M\^) = sin^9. Such a metric tensor 
is also what is required in order to compute the scalar quantity we call the norm of 
a (tangent) vector, e.g., the norm of the velocity vector dx/dt is given by 



Similarly, one can compute the angle between two (tangent) vectors, e.g., two veloc- 
ities, at each position in the x-space. 

The situation in the configuration space of a fluid with coordinates (p, T, u) is totally 
different from that in the x-spaces described above as long as it is not known how 
to introduce, in a natural and intrinsic way, a geometric structure, e.g., something 
similar to the metric tensor M\^{x) in (|l]) and @. Such a lack of geometric structures 
imposes many restrictions on the kind of computations that one can perform: it is 
impossible to talk of "the distance" between two states of the fluid, i.e., between two 
positions in the space (p, T, u) ; neither is it possible to talk of "the norm" of the 
vector formed by the rate of change of the dynamical variables {dp/dt, dT/dt, du/dt); 
it is not possible to consider the angle between two such vectors; there is no volume 
element defined in configuration space, therefore, it does not make sense to talk about 
"the density" of a distribution of points in that space; etc. These limitations are too 
restrictive since, to name but a few examples: 1) A norm is needed when studying the 
(in)stability of flows, especially when the linear growth of the perturbations is only 
transient 0; 2) The angle between two directions is required in order to determine 
how strongly a given perturbation projects along an optimal perturbation [Q; 3) A 
volume element is needed in order to determine the density distribution of ensemble 
simulations of flows Therefore, it is often opted to circumvent these limitations 
by introducing an acceptable, if somewhat arbitrary, metric tensor, e.g., the kinetic- 
energy metric in the case of incompressible, isentropic flows 0. Such a situation is 
not unique to fluid mechanics: often our knowledge of a dynamical system is limited 
to the variables that define its phase space and to the equations that determine the 
time evolution of these variables while any information on the geometry of the phase 
space is lacking. 

In this article it is shown how to generate, in a natural way, not only the ana- 
logue of the spatial metric tensor Mxi,{x) but also a rich geometric structure in the 
configuration space of an inviscid, moving fluid starting from the fluid's probability 
distribution. The derivation of these results makes use of ideas and techniques from 
statistical mechanics, see e.g., P, and from statistical manifolds theory, see e.g., 
0, a particular branch of differential geometry, see e.g., 0, ^. The presentation 




dx dx' 



.V 



(4) 
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is as self-contained as possible. A complete description of the required background 
material, however, was not attempted; detailed bibliographic references are given for 
the benefit of those readers willing to go deeper into some aspects. 
The paper is structured as follows: In Section 2, the probability density of a moving 
fluid in local thermal equilibrium is presented. In Section 3, some fundamental con- 
cepts of the theory of statistical manifolds are introduced. In Section 4 these ideas 
are applied to a moving fluid, the natural metric tensor in the fluid's configuration 
space is derived and an application is described; some additional information is given 
in Appendix A. In Section 5, it is shown that the probability density generates a fiat, 
but non-Riemannian, geometry in the configuration space. Basic concepts on flatness 
(of affine connexions) are reviewed in Appendix B. Based on this flatness, a distance- 
like function is defined in Section 6. The computed metric and distance-like function 
are Galilean invariant, as they should. When the thermodynamic variables can be ig- 
nored, both the metric and the distance-like function reduce to the above-mentioned 
kinetic-energy norm times a conformal factor. For the purpose of illustration, the 
expressions corresponding to an ideal gas are explicitly given in Sections 4 and 6. In 
Section 7, we comment on related work done by others; in Appendix C, the differences 
between the approach developed in the present paper and a particular Hamiltonian 
formulation of fluid dynamics are pointed out. Finally, in Section 8, we review the 
results, discuss their generalizations, e.g., to systems far from thermal equilibrium. 



2 Statistical mechanics of a moving fluid 

The first step in our construction of an intrinsic geometric structure in the (p, T, u) 
space consists in associating a probability density to these variables. In this Section 
we derive the fluid's probability distribution by applying statistical mechanics to a 
moving fluid in local thermal equilibrium. 

Consider a simple fluid characterized, on a macroscopic level, by its density field 
p{x,t), temperature field T{x,t) and velocity field u{x,t); from here onwards, we do 
not express the (x, t)-dependence explicitly. We adopt the standard assumptions that 
make possible the derivation of the Euler equations and, less rigorously, of the Navier- 



Stokes equations from local thermal equilibrium statistical mechanics, see, e.g., [ITU 
Chapters 1 to 3 of Part 1] and |[Tl|, Chapters 1 to 3], to wit: One considers systems in 
which the above-mentioned fields vary on length scales that are much larger than the 
intermolecular distances and on time scales that are much larger than the microscopic 
time scales needed for local thermal equilibration. Then the space can be divided into 



3 



small volume elements AV with a typical size much smaller than the length scales 
of the macroscopic fields but much larger than the intermolecular distances so that 
the statistical mechanical description applies to them, i.e., on the length scales of 
A^, the fiuid is in local thermal equilibrium. The dynamics of the fiuid is assumed 
to be such that the following extensive quantities are conserved: a) the number of 
(indistinguishable) particles A^, b) the total momentum of these particles 

N 

M:=mJ2vk, (5) 

k 

where m is the mass of the particles, Vk is the velocity of the A;-th particle, and c) 
their total energy 




(6) 

Vki being the potential interaction between particles k and /. (The symbol to the left 
of := is defined by the expression to the right.) 

Let us denote the conserved quantities by {ifj(^Ar) \i = 1, . . . , s} where := {xi, . . . , 
xjy, mvi, . . . , mviy \N = 1, . . . , oo} and the values of {Hi} in a thermal equilibrium 
state by {rji}. These two elements, i.e., the expressions defining the conserved quan- 
tities {Hi} as functions of the microscopic variables {^n} and the values {rji} charac- 
terizing the macroscopic state, are the essential ingredients needed for the statistical 
mechanical description of the system. Statistical mechanics establishes that the prob- 
ability density in the phase space C,n, of a system in AV, in thermal equilibrium with 
and free to exchange particles, momentum and energy with a surrounding thermal 
bath is 

pi^M-e) := exp [e^H,{^N)]/c^mz{e), (7) 

see, e.g., Chapter 4], Chapters 9, Section B.3]. We use the standard short- 
hand conventions: repeated lower and upper indices are summed up; 6 stands for 
{6^\i = 1, . . . ,s}. The 2{6) in the denominator, called the grand partition function, 
is the normalizing factor needed to make sure that J2n IdC,Np{C,N', G) = 1, i-e., 

oo „ 

Z{9) := ^ / exp [e'H,{^N)]/c'' N\, (8) 

the constant c makes each of the terms contributing to Z dimensionless. From (0) 
and (^ we see that the average values of the conserved quantities {ifj(^Ar)}, which 
we call {rji}, can be written as 

d\nZie) 

V^{0)= g^y , z = l,...,s. (9) 
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Since the thermal equihbrium state of a system with s conserved quantities is com- 
pletely characterized by their macroscopic values {rji}, the last expression makes 
evident that the whole thermodynamics can be obtained from the grand partition 
function Ziff). In this way, one identifies In^(^) as the thermodynamic potential 
which is minimized at fixed values of the 6'- variables and the 6'-variables as the inten- 
sive thermodynamic parameters (like temperature, pressure, chemical potential, etc). 
Since the system in !W is in thermal equilibrium with its surroundings, the values 
of these intensive variables must be equal to those of its thermal bath, see, e.g., 
Chapter 4], []^, Chapters 9, Section B.3]. 

In the case of a moving fluid, the conserved quantities are {N,E,M), confer and 
(I), so that s = d + 2 and the probability density (|^ reads 

p{N, E, M; 7, -(3, k) = exp[-fN -(3E + K- M]/c^N\Z{-f, -f3, k), (10) 

i.e., 9 = (7, — /?, k). It follows then that the average, macroscopic values of (A^, E, M) 
are given by 

(AT) = ^^"^^2'"^'"^ (11) 

(77 

(E) = _^ili£|LzM) (12) 
(M.) = ^'°y-^^' . A=l.....d. (13) 



where the pointed brackets indicate average over the probability distribution (|iy). It 
is convenient to introduce the macroscopic velocity of the fluid u by 

m{N)u := (M) (14) 

and the specific internal energy 6(7, /3, k) as 

{N){e + -u'):={E). (15) 

Relating these expressions to the thermodynamics of the system, one identifies In Z{'j, —f3, k) 
as the thermodynamic grand potential, i.e., 

\nZ{-f,-(3,K) = I3AVP, (16) 

where P = F(7, P, k.) is the pressure, P as the inverse local temperature, P = {kT)~^, 
k is Boltzmann's constant and 7 as 7 = with /i the local chemical potential; see, 
e.g., [^, Chapter 4], Chapters 9, Section B.3]. Due to the simple dependence of M 
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and of E upon the particles' velocities Vk, confer (j^) and the integrals over the 
velocities in (|]) can be performed; from (0) and ([T^ one obtains then that 



K = pu (17) 

and that 



Z(7,-/?,«:) = Z(7-^,-A0). (18) 



Galilean invariance implies then that 



7 = 7 + ^, (19) 



where 7 is the value of 7 for the same system at rest. Making use of (|Ty|), one sees 
then that /i = fl + mu'^ /2, i.e., the chemical potential of the fluid moving with velocity 
u is shifted by an amount mv? /2 with respect to /i, the chemical potential of the fluid 
at rest. 

In the theory of statistical manifolds 0] it is shown that a probability density like 
p{^N]d) in (0) induces a geometric structure in the parameter-space 9, i.e., in the 
space of the intensive thermodynamic variables. This is worked out in the following 
Sections. 



3 Statistical manifolds 

Given two probability densities, say p{C,N, Oi) and p{C,N, O2), how should one quantify 
their (dis)similarity? Or given three distributions corresponding to 61,62 and ^3 
respectively, is it possible to determine which pair of distributions is "closer" than 
the two other pairs? These and related questions have been extensively studied in 



statistics, see, e.g., [|T^, pages 290-296] and ||T3[. In this Section we give a fleeting 
overview of some concepts and results from the theory of statistical manifolds that 
we apply later to a moving fluid. 

If 0(61,62) is a scalar quantity measuring the difference between two probability 
densities p{^N, 61) and p{C,n, 62), then it should satisfy: 

1) D{6i, 62) > 0, the equality holding if, and only if, p{^N, 61) = p{iN, 62), 

2) -0(6*1, 62) is sufficiently continuous in p{^N, 61) and p^^n, 62)- 

In the statistical literature such measures are often called divergences or contrast 
functionals; in order to avoid any possible confusion with the usual mathematical 
meaning of divergence, we shall keep to the second name. 
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We shall impose two additional constraints: 

3) The value of D{9i,92) should be independent of our choice of coordinates in ^at- 
space, i.e., it should be invariant under general coordinate transformations in the 
space of the random variables. 

4) We shall consider only contrast functionals that are local in the random variables 
{^n}^ 

In order to satisfy these constraints, it is sufficient that the D{6i,62)s be of the 
following form 0, [0, pp. 349-350], 



BW.«,) = ^Up&,«.)/(^), (20) 

where the function / must be sufficiently smooth and should satisfy /(I) = and 
/(t) — — 1) > if t 7^ 1. Without any loss of generality, one normalizes / such 

that /"(I) = 1. 



Consider next two distributions that are infinitesimally close, i.e., 6 and 6 + dO. It 
turns out that, up to third order in dO, the corresponding contrast is given by [|I], 
theorem 3.10] 

D{e, e + de) = ^ g,,{e) dO' de^ + ^ j; k] + |T,,fc(e)| de' dO' de'' + oide^), (21) 



where the components of gij{0), known as the Fisher tensor |T^, |T3[, are given by 



the pointed brackets indicate an average taken over the p{Cn,0) distribution. From 
this expression it is evident that gij is symmetric in the indices i and j and that, 
since it is an addition of products of covariant vectors, under a general coordinates 
transformation it transforms like the product of two covariant vectors, i.e., it is a 
covariant tensor^. Moreover, any sensible choice of the parameters 6 ensures that it 
is non-singular. 

The first term in the coefficient of the third-order contribution to the expansion pT|), 
[i,j;k], is given by 

U rk]--(^ + ^- ^] (23) 



^C.R. Rao has considered also non-local contrast functionals, ^, pp. 226-227]. 
^For a more detailed meaning of these geometric concepts, see, e.g., [|[ ^. 
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its meaning is discussed in Appendix B and in Section 5. The second term in the third- 
order coefficient is a totally symmetric covariant tensor, often called the skewness, 

Both symmetric tensors, gij and Tjjfc, will play an important role in the sequel. 
The constant a in (pT]) is given by a := 2/'"(l) + 3. 

A number of points are worth noticing: 

• The symmetric contravariant tensor gij{6) has all the credentials for being the 
natural metric tensor in the ^-space. Besides the clear meaning ensuing from (plj), it 
has an important statistical significance |]15|, |16] that is described in Appendix A. 

• Up to second order, the expansion ( p]) is independent of the particular choice of 
the function / that defines the contrast D, confer (PO]). (Remember that / has been 
normalized so that /"(I) = 1.) 

• Up to third order, the whole dependence upon the function / has been brought 
back to a single number, namely, to a. 



4 Fisher's metric 

In this and the following Sections, we apply the theory of statistical manifolds to the 
probability density p{C,n, 0) of a moving fluid. 

As shown in the previous Section, the Fisher tensor appears naturally as the metric in 
^-space. From (|TOD it follows that in the coordinate system (6'^, . . . , 6'''^'^) = (7, — /3, k), 
the metric tensor reads, 

dHnZ{9) ,,,d^(3P 



see also (^) in Appendix A. We see then that the metric elements are expressed in 
terms of the derivatives of the thermodynamic potential. Moreover, from (pl^-[T3|) or 
from (^) one has that, in the (7, — /?, R) coordinate system, 

dr]i = gij{9)de\ (26) 

where rji are the average values of the extensive quantities given in (|lT]-0). I.e., in 
this coordinate system, the lowering of the indices corresponds to passing from the 
intensive variables 6 to the extensive ones 77. The inverse transformation is achieved 
by raising the indices by means of the inverse metric g^^{9). 
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It is instructive to work out the metric tensor for the case of an ideal gas, i.e., for a 
vanishing intermolecular potential Vki in (P). In this case all the integrals over the 
molecular velocities Vk and over their positions in can be performed and one 
finds that 

lnZ(7,-A/?) = Al^r'/'exp(^-7). (27) 

It is convenient to express the results in the coordinate system {p,f3,u), where p is 
the mass density, i.e., p := m{N) / AV . The non-vanishing elements of the metric are, 

9pp = AV{mp)~^, 

gpp = AV-—r-, 
2 mp^ 

gu^u, = AVp(35,,. (28) 

Therefore, each volume element AV contributes a squared distance {dV)"^ between 
two states (p, /5, u) and (p + dp, [3 + d(3, u + du) given by 



AVp 



m 



— ) +f (-tI +m[3du-du 



(29) 



This expression agrees with our expectations: 1) The metric coefficients are inde- 
pendent of u, as demanded by Galilean invariance^; 2) Since no prefered or external 
scales are available and taking into account the Galilean invariance, the differentials 
must always appear in the form of dp/p and df3 / (3\ the gas temperature determines 
the scale for measuring the kinetic energy associated with du\ 3) Cross-terms like 
dp du and df3 du cannot appear due to the rotational symmetry of the gas and 4) The 
factor multiplying the square brackets is just the number of particles in AV . 
What could not have been guessed beforehand, is the absence of the cross-term dp d(3 
and the precise form of the coefficients. 

One of the simplest applications of a metric tensor is the computation of the norm 
of a tangent vector, as in (p. Now we can do this in the configuration space of the 
moving fluid, i.e., we can compute the norm of the local rate of evolution of our sytem, 
that we denote by F, as 

:= g^jO'e^, (30) 

where the dots indicate the time derivatives of the corresponding quantities as given 
by the Euler equations. In order to distinguish between this generalized velocity F 
and the standard velocities Vk and u, we shall call this quantity the rapidity. The 



^The Galilean invariance of the metric in the general case (non-ideal gas) is proven in Section 6, 
confer (46). 
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rapidity squared is a scalar with dimension [time] ^. From ( ]26|) we see that it can 
also be written as 

= m9\ (31) 

i.e., it is the contraction between the velocities of the intensive variables with the 
velocities of their corresponding conjugate extensive variables. 

Let us compute the rapidity associated with the material derivatives in the Euler 
equations with no external forcing, i.e., 

dp ^ 
— = -pdivM, 
at 

pCv—J- = — -T dlVM, 



dt 



du - „ 



(32) 



where is the specific heat at constant volume. In the case of an ideal gas, i.e.. 



making use of (|28|) and inserting (|32D into (|30|) leads to 



m 



1 + ^ 1 (divM)2 + B!T Vln(pT) ' 



(33) 



where R' := k/m is the specific constant of the gas. 



5 Flatness of the manifold 

Another standard application of a metric is the determination not only of an infinites- 
imal distance, as in (pOj), but also of finite distances, say L(l,2) between positions 
1 := (pi, — Ml) and 2 := (p2, —(32, U2). This distance is defined as the length of the 
shortest path connecting 1 with 2, 



Lil,2)=mmj'dT^g,,i9)^^, (34) 

where the minimum is taken over all paths from 1 to 2 and r is any parametrization 
of these paths. 

This calculation is not a simple one when there is no coordinate system in which the 
metric tensor gij{6) takes a simple form; this happens in the case of our moving fluid 
since, as we will see in this Section, (|55| ) and comments thereafter, (^) is not fiat. 
Luckily, it turns out that the geometry of statistical manifolds offers more interesting 



and useful possibilities than (|34 ); this is particularly so in the case of a probability 
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density of the form (0) as is our case. In order to reahze this, we need some concepts 
from differential geometry; for the sake of completeness, these concepts are listed in 
Appendix B. 

Two simple calculations indicate that the connexions described in Appendix B may 
play an important role in the case of our moving fluid: 

Let us compute first the [ij;k] connexion compatible with the Fisher metric p5|); 
from (|50D one has that, in the coordinate system {6^, . . . , 6'^'^'^) = (7, —j3, k), 

2 06^06^6'^ ■ ^^^^ 
One can compute the corresponding curvature tensor and check that, even for the 

ideal gas metric (|28|)-(p9D, this connexion is not flat. 

Next, let us compute the skewness tensor Tij^ generated by (|TD|), confer (P^); one 
finds that 

^'^'^^^ - Wdew- ^^^^ 

From these two last results, we see that by taking [ij; k] — {l/2)Tijk one gets a con- 
nexion that vanishes identically in this coordinate system, i.e., this connexion is flat. 
It should be noted that the existence of a fiat connexion implies that there is a pref- 
ered family of coordinate systems, in our case, the system of the intensive variables 
9 and their linear combinations. The physical relevance of these variables becomes 
even more evident when one considers the phenomenon of phase-coexistence, e.g., 
liquid-vapour coexistence: coexisting phases correspond to fiat portions of the ther- 
modynamic potential surface only in terms of the intensive variables. 

Similarly, one can check that also [ij\ A;] + (l/2)Tjjfc, notice the change in sign, vanishes 
identically in the extensive variables coordinate system 77, i.e., also this connexion is 
fiat. Connexions of the form [ij; k] — {a/2)Tijk play an important role in the theory 
of statistical manifolds; they were introduced by Chentsov fl^, Efron |18| and Amari 
Chapter 3]. In the third reference it is shown that if such a connexion is flat for 
a certain a then it is also fiat for —a. In fact, the coordinate systems 6 = (7, —f3, k) 
and 7] = {{N), (H), (M)) play totally symmetric roles: it can be shown ^ Theorems 
3.4 and 3.5] that the inverse of the metric tensor, is also given by the second-order 
derivatives of a function ^{rj), this time with respect to the extensive variables rj, i.e., 

- a^' (^^) 

that the function $(77) is nothing else but the Legendre transform of In Z{6) and that 
the coordinate systems are related to each other as in a Legendre transformation, i.e., 

(9$ 

0\V) = T^, (38) 
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^(^) + $(r7) -6'-r]i = 0, 



(39) 



where we have introduced \1/ := In 2, confer also (P). Recall that all functions are 
assumed to be sufficiently smooth so that ( p9| ) is good enough for the definition of 
the Legendre transform; for more complicated situations, see . 

We have shown then that the probability density of the moving fiuid generates not 
only a metric tensor (|25| ) but also two fiat connexions, namely [ij;k] ± {l/2)Tijk 
and that all this is closely related to a Legendre transformation of thermodynamic 
potentials and variables. Some interesting consequences of these facts are presented 
in the next Section. 



6 Distance-like function 

While the computation of the distance L in (0) implies difficult calculations, it was 
pointed out by Amari Section 3.5] that, when the connexions [ij;k] ± {l/2)Tijk 
are fiat, it is possible to define a distance-like function D{1,2) between positions 
1 := (71, Ki) and 2 := (72, -/?2, '^2) as follows 

D{1, 2) := + $(2) - 9\1) ■ r],{2), (40) 

with ^ := InZ as in (||), confer also and (|). 

Amari has shown that this function shares some essential properties with the usual 
distance functions, to wit"^: 



1) 1^(1, 2) >0, 
the equality holds when, and only when, 1 = 2, (41) 

2) D{9,e + de) = 
^ gij{e) dO' de^ + ^ j; k] + ^T,,fc(^)| dO' dO' d9' + 0{d9'), (42) 

3) Dil,2) = D{l,Q) + D{Q,2), 
where Q is connected to 1 by a 0— geodesic and to 2 by an 
77— geodesic and these two geodesies intersect orthogonally at Q (43) 

and 4) The mmD(l,Q) 
is obtained at a point on the boundary of the smooth closed region Q 



*The first property is a special instance of Fenchel's theorem [n9| 
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that is the projection of 1 along a 6'— geodesic orthogonal 

to the boundary of Q. (44) 

The ^-geodesies above is a linear interpolation from ^^(1) to 6{Q) while the r^-geodesic 
is a linear interpolation from r]{Q) to 1]{2). The third property is a generalization 
of Pythagoras' theorem to a space with two biorthogonal coordinate bases, related 
to 6 and rj in our case. Analogously, the fourth property generalizes the notion of 
projection to such a space. The reader should refer to [0, Section 3.5] for the proofs 
of the properties listed above. 

There is one important difference between the usual distance functions and the D 
defined above: In general, the D function is not symmetric, i.e., D{1, 2) ^ -0(2, 1). 
This asymmetry is associated with the fact that the path going from 1 first along an 
77-geodesic to Q' and then along an orthogonal ^-geodesic to 2 is, in general, different 
from the path described in (|43|). If one is set on defining a symmetric distance, 
then taking the minimum of the values D{1,2) and D{2, 1) seems to be the most 
satisfactory solution. 

From properties ( Pl) and (^21) , we recognize that this D may belong to the family 
of contrast functionals discussed in Section 3. In fact, using (^, (p5D, ( |38D and (p9|), 
one finds that D{1,2) in (^Oj) can be written as 

D{1,2) = {[32 -f3i){H)i- {k2 - ■ (M)i - (72 - 71) {N)i + AV{f32P2 - APi), (45) 

where indicates that an average is taken over p{N,E,M;6i). One can check 

then that, in our case, the last expression for -0(1, 2) is identical to ( ^0]) with f{z) = 
z In z. The circle is now complete. 

The Galilean invariance of D{1,2) follows from introducing ([T5|), (|17|) and (|l^) into 
the last expression above; one finds then that the contrast functional (^) or (^) can 
be written as 



m 



D{1, 2) = {(32 - (3i){N),e, + -(32\u2 - Uir{N)i - (72 - 7i)(iV)i + AV{(32P2 - (3iP, 



(46) 

Combining this with ( pTl) and with (^2|) , one sees that both gij{6) and Tijk{6) are 



always independent of u, i.e., Galilean invariant. 

Notice also that when the variations in the thermodynamic variables can be ignored, 
as is the case in an (effectively) incompressible, isentropic flow, the contrast reduces 
to -0(1, 2) = {AV pP/2)\ui — U2\'^ with p = pi = P2 and /? = /?i = P2, i-e., the kinetic- 
energy norm times a conformal factor. 
In the case of an ideal gas, the contrast function reads: 



0)(1,2) = —AV 

m 



—P2U2-U1 + ln— + l +- ln— + — -1 

2 V P2 Pi / 2 \ ^2 A , 



(47) 
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7 Related work 



Before closing, we review and comment some papers that deal with related questions: 
Weinhold proposed to use the matrix of second-order derivatives of the internal 
energy as a metric in the space tangent to the equation-of-state surface at an equilib- 
rium point. Some researchers tried then to attach a physical meaning to the curvature 
tensor derived from the Weinhold's metric-compatible connexion, confer Section 5. 



This led to a lengthy discussion which has been summarized in all these re- 
searchers ignored the fiat connexions [ij; k] ± |Tjjfc and the associated distance-like 
function D as discussed in Section 6. In the context of the present article, Weinhold's 
metric comes close to the Fisher metric (^) while Gilmore's approach, see is 
closer to the distance-like function D (|40|) . 

The Fisher metric was introduced into thermodynamics by Ingarden |2^, Janyszek 



and Mrugala p4| studied the curvature tensors of the corresponding metric-compatible 
connexions and tried to associate this curvature with some physical properties; they 
did not consider the fiat connexions [ij; k] ± ^Tijk neither the associated distance-like 
function D. 

On the 18th of December 1995, I presented this paper at the Technical University 
of Berlin. Prof. U. Simon informed me then about affine and projective differential 
geometry which deals, among other things, with conjugate connexions like [ij; k] ± 
|Tjjfc and their flatness. Excellent, up-to-date reviews of this branch of differential 
geometry are [25| and 



There exists a rich literature on the Hamiltonian (symplectic) structure of hydrody- 
namics that can also be seen as one branch of differential geometry; see, e.g., the 
review articles the references therein and p9| , ^ for another interesting ap- 

plication. This approach is based on the Poisson bracket and it does not require a 
metric tensor. By contraposition, the present article does not use the Poisson bracket 
and leads to a metric tensor, flat connexions, etc; i.e., the two approaches are not 
necessarily related and can be seen as complementary. In Arnold's analysis of incom- 



pressible flows ||31|, a metric plays an important role, however, this metric is just the 
metric of the ambient d-dimensional space that appears in the kinetic energy term, 
i.e., is not the type of metric developed in the present article. There is one instance of 
a Hamiltonian structure for fluid dynamics that apparently contains a candidate for 
a metric tensor in the configuration space of a moving fluid; this instance is discussed 
in Appendix C. 
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8 Summary and discussion 



In this article, it has been shown that tlie grand canonical partition function (p^) 
which describes the thermal fluctuations of a moving fluid in local thermal equilibrium 
generates a natural metric tensor QijiO), two flat connexions [ii',k] ± \Tijk and the 
corresponding distance-like contrast function _D, confer (p5|), (^), (p6D and (^0|). In 
the case of an ideal gas, the explicit expressions for these quantities have been given 
in ( p9D and in (^T]). It was also shown that these quantities are Galilean invariant, 
as they should, and that they reduce to the kinetic-energy norm times a conformal 
factor when one can ignore variations in variables other than the velocity field, confer 
(^). Noteworthy aspects of the geometric structure are that the flat connexions 
are not metric-compatible and that the distance-like function D is not symmetric. As 
pointed out by one of the referees, the analysis presented in this article can be applied 
to relativistic fluids, to plasmas, etc. In fact, some of the results listed above, e.g., the 
particular form (|25D of the metric tensor, as well as (|37D -(p9D and the flatness of the 
connexions [ij; k] ± ^Tj^fc, are valid for all systems characterized by Gibbs' probability 
distributions like (0). 

The following considerations are relevant with regard to further applications of our 
approach. We have seen that all the components of the geometric structure are gen- 
erated from \nZ = AV(3P and since, in the absence of external forces, the equation 
of state P = P(7, (3) is all that is needed in order to write down the Euler equations 
of fluid motion, it is then clear that the metric, connexions and contrast function are 
closely related to the dynamical equations. When the system is far from thermal equi- 
librium, however, a different approach may be more appropriate. For example, when 
considering a fluid in a turbulent stationary state, the corresponding probability den- 
sity should be used. Similarly, when studying the predictability properties of systems 
in the presence of noise, a time-dependent probability density should be employed 
132| . In this sense, the geometric structures we have described are not universal; the 



specific physical conditions of the system under consideration determine whether a 
thermal-equilibrium, a stationary-state or a time-dependent probability density is the 
most appropriate starting point. 
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Appendix A: On Fisher's metric 



Besides the clear meaning due to (pi]), the Fisher metric tensor (^) plays an impor- 



tant role in statistics due to the following theorem |p!5| , |T6|: 

— * 

Suppose that we perform n independent measurements of the random variables (A^, E, M) 
and use them in order to estimate the value of the parameters (7, — /?, k) in ([T0|) . De- 
note by {&'■} an unbiased estimate of these parameters. Then, for large enough n, the 
covariance of these estimated values with respect to the exact values {6^} = (7, — /3, k) 
has the following lower bound 



cov 



{&' - e'){e^ - e^)] > -g'\ (48) 

-I n 



where {g^^) is the inverse of the matrix {gij). In other words, given a fixed tolerance 
error, the larger the distance between two probability distributions as measured by 
gij, the smaller the number of measurements needed in order to distinguish between 
them. 

Finally, one should notice the following identity, 

, , ld\n.pd\n.p\ / d^\np\ , , 

This identity is obtained as follows: Jdxp{x,6) = 1 implies that {dlnp{x,6)/d6^) = 
Jdxdp{x, 6)/d6^ = 0; taking then the derivative of {dlnp{x, 6)/d6^) with respect to 6^ 
leads to the identity in (^91). (It is assumed that the order of derivation with respect 
to the {6'*} and integration over x can be interchanged.) 



Appendix B: On afRne connexions 



In some cases, the geometry of a manifold is not completely defined only by its metric 
tensor but also by another geometric object, called the connexion of the manifold^, 
see, e.g., [§, Chapter 4, Sections 28 and 29] and Chapters I-IX]. Briefly stated, the 
connexion expresses mathematically what it means "to move a vector along a curve 
in such a way that the vector remains constant", i.e., it defines the parallel transport 
of a vector. 

We list some important properties of the connexions that we need in Section 5: 

• Under a general transformation of the coordinates, a connexion does not transform 

like a tensor, i.e., a connexion is not a tensor. In particular: a symmetric connexion 

^Whcn this happens, we speak of a non-Ricmannian manifold. 
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may vanish identically in a particular coordinate system without vanishing in a differ- 
ent one while a tensor that vanishes identically in a particular system of coordinates 
does so in all coordinate systems. 

• The result of adding a third-rank tensor to a connexion is a connexion. 

• We say that a manifold is flat with respect to a connexion when every vector 
that is parallel transported along every closed path returns to its original condition; 
otherwise, one says that the manifold is curved. 

• A manifold is flat with respect to a connexion if, and only if, it is possible to find 
a coordinate system in which this connexion vanishes identically. 

• Whether a given connexion is fiat or not can be determined in a coordinate-system 
independent way: a connexion is fiat if, and only if, it is symmetric in its two first 
indices, i.e., torsionless, and its curvature tensor vanishes identically. 

• The partial derivatives of a tensor are not tensors, however, using the connexion 
it is possible to define the so-called covariant derivatives which are tensors. If the 
connexion vanishes identically in a given coordinate system then, in that coordinate 
system, the covariant derivatives coincide with the partial derivatives. 

• There is only one connexion symmetric in its two first indices such that the co- 
variant derivatives of the metric tensor vanish identically. If Qij denotes the metric 
tensor, then this connexion^ is given by 

1 dgjk _ dgij\ 

2 \ dei 89^ 86'- ) ■ ^ ' 

This is precisely the [ij ; k] appearing in the third-order term in PT| ) , confer (^) . This 
connexion is called the metric-compatible or Levi-Civita connexion. For the proofs 
of these statements, refer to, e.g., P, Chapter 4, Sections 28 and 29] and Chapters 
I-IX]. 

Appendix C: On the Dubrovin-Novikov approach to flows 



In an illuminating article , Dubrovin and Novikov introduced a novel Hamiltonian 
formalism for one- dimensional systems of hydrodynamic type described by s dynam- 
ical fields. (An extensive review of related ideas and developments can be found in 
P^ .) One notable aspect of this formalism is the central role played by a symmetric 
covariant tensor of type (2,0), call it A*-', 1 < (i, j) < s: when det(A*-') 7^ 0, the 
Jacobi identity that the Poisson bracket has to satisfy implies that the Levi-Civita 
connexion generated by Ajj := (A^^)jj must be fiat. Moreover, when det(A*-') 7^ 
and other, relatively mild conditions are satisfied, it is possible to reconstruct (up 

^Actually, this connexion is obtained by raising the A;-index of the [ij; k] in (pO[), i.e., by g^'lij; k]. 



17 



to a constant factor) the tensor A*-' directly from the hydrodynamical equations of 
motion [^. Therefore, in such cases, it is possible to associate with the equations of 
motion, one metric in the phase space of the fluid: the inverse of (A*-'). In practice, 
these symmetric tensors cannot be used as metrics but in few, unrealistic cases due to 
the following reasons: 1) When the fluid exists in a (i-dimensional space, d different 
symmetric tensors A*-' must be introduced [P3] leading, in the best case, to d different 
metrics in the fluid phase space; 2) Even in the case of a one-dimensional physical 
space, one often has that det(A*-') = 0, this is the case, e.g., of a one- dimensional 
nonbarotropic fluid PSI p. 59]. Moreover, it should be noticed that the A*-' 



tensor is part of the Poisson bracket definition, therefore, different dynamical systems 
generated by different Hamiltonians, i.e., fluids with different equations of state, but 
by the same Poisson bracket share the same A*-''; this property is obviously not de- 
sirable when one is trying to identify the geometric structures that characterize and 
distinguish between different hydrodynamical systems. Similarly, some hydrodynam- 
ical equations have a multi-Hamiltonian formulation, i.e., the same set of equations 
can be generated from different Hamiltonians by means of different Poisson brackets 



37, 3Sl; in such cases, each of these Poisson brackets would lead to a different metric 



tensor for the same dynamical equations. 
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